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ABSTRACT 

We present early- and late-time photometric and spectroscopic observations 
of supernova (SN) 2009hd in the nearby spiral galaxy NGC 3627 (M66). This SN 
is one of the closest to us in recent years and provides an uncommon opportunity 
to observe and study the nature of supernovae. However, the object was heavily 
obscured by dust, rendering it unusually faint in the optical, given its proximity. 
We find that the observed properties of SN 2009hd support its classification as a 
possible Type II-Linear SN (SN II-L), a relatively rare subclass of core-collapse 
supernovae. High-precision relative astrometry has been employed to attempt to 
identify a SN progenitor candidate, based on a pixel-by-pixel comparison between 
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Hubble Space Telescope (HST) F555W and F814W images of the SN site prior 
to explosion and at late times. A progenitor candidate is identified in the F814W 
images only; this object is undetected in F555W. Significant uncertainty exists 
in the astrometry, such that we cannot definitively identify this object as the SN 
progenitor. Via insertion of artificial stars into the pre-SN HST images, we are 
able to constrain the progenitor's properties to those of a possible supergiant, with 
intrinsic absolute magnitude Mp 555W ^ —7.6 mag and intrinsic color (V — I) ^ 
0.99 mag. The magnitude and color limits are consistent with a luminous red 
supergiant; however, they also allow for the possibility that the star could have 
been more yellow than red. From a comparison with theoretical massive-star 
evolutionary tracks which include rotation and pulsationally enhanced mass loss, 
we can place a conservative upper limit on the initial mass for the progenitor of 
M ini <^ 20 M Q . If the actual mass of the progenitor is near the upper range allowed 
by our derived mass limit, then it would be consistent with that for the identified 
progenitors of the SN II-L 2009kr and the high-luminosity SN Il-Plateau (II-P) 
2008cn. The progenitors of these three SNe may possibly bridge the gap between 
lower-mass red supergiants that explode as SNe II-P and luminous blue variables, 
or more extreme red supergiants, from which the more exotic SNe Il-narrow may 
arise. Very late-time imaging of the SN 2009hd site may provide us with more 
clues regarding the true nature of its progenitor. 

Subject headings: galaxies: individual (NGC 3627) — stars: evolution — super- 
novae: general — supernovae: individual (SN 2009hd) 



1. Introduction 



The Type II-Linear supernovae (SNe II-L) are among the least common and, there- 
fore, most poorly studied subclasses of core - collapse supern ovae (CC-SNe). They represent 
6.4-10% of all CC-SNe JSmith et al l l201l[ iLi et al.l lioilh . The spectra of SNe II-L are 
similar to those of the much more common SNe Il-Plateau (II-P), but SNe II-L are dis- 
tinguished by the shapes of their light curves. While light curves of SNe II-P exhibit a 
relatively constant plateau for ~100 d after an initial peak in the first few days, the light 
curves of SNe II-L instead show a linear decline commencing just after the initial peak (e.g., 
Barbon. Ciatti. fc Rosinolll979f ). The very short, or nonexistent, plateau phase for SNe II-L 
may well result from a low- mass H envelope at the time of explosion (between 1 and 5 M ; 
Swartz. Wheeler, fc Harknesslll99ll ). The envelope can be further reduced through pre-SN 
mass loss, and if so, we would infer that the mass of the progenitor star is higher than that 
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expected for SNe II-P, 8.5lJ;g £ M ini £ 16.5 ± 1.5 flSmartt et al .1120091) . The direct 
identification of the progenitor of the SN II-L 200 9kr (|Elias-Rosa et al. 2010 ) provides evi- 
dence that this indeed may be the case (note that Fraser et al. 20091 suggest that this SN is 
a spectrally peculiar SN II-P). Clearly, the identification of the progenitors of other SNe II-L 
is necessary to establish a trend, if one exists. 

In this paper, we present the case of SN 2009hd (a = ll h 20 m 16 i ;90 S = +12°58 / 4 7"1; 



J2000.0) in NGC 3627 (Messier 66, M66). This SN was discovered bv iMonardl (120091 ) on 
2009 July 2.69 (UT dates are used throughout) and was spectroscop ically classified as a 
reddened SN II (IKasliwal. Sahu. fc Anupamall2009l) with slow ev o lution ( Kasliwalll2009l) . and, 
subsequently, as a SN II-P (IBerger. Foley. &: Covarrubiasll2009l ). iMaoa (120091 ) first identified 
a possible progenitor star in high spatial resolution archival Hubble Space Telescope (HST) 
images obtained between 1997 and 1998. Here we show that both the photometric and 
spectroscopic observations suggest that this object is, in fact, a SN II-L. We also present a 
substantially improved identification of the progenitor, based on a pixel-by-pixel comparison 
between the pre-SN HST data and late-time HST ima ges of the SN. We not e that M66 has 
been host to t wo other known S Ne, the SN II 1973R (ICiatti fc Rosindll9770 and the Type 



la SN 1989B (IWells et all 119941 ) . as well as the "SN impostor" SN 1997bs (IVan Dvk et al. 
2000l ): see Figured) 



2. The Nature of SN 2009hd 
2.1. Photometry 

Optical (BVRI) images of SN 2009hd were obtained with the 1.3 m SMARTS+ANDICAM 
at Cerro Tololo Inter-American Observatory (CTIO), the Palomar 1.5 m telescope, the 2.0 m 
Liverpool Telescope and the 2.6 m Nordic Optical Telescope (NOT) at the Roque de los 
Muchachos Observatory, the 6.5 m Baade Magellan telescope, and the 10 m Keck-II tele- 
scope. The photometric observations were processed using IRAF0 routines with the stan- 
dard recipe for CCD images (trimming, overscan, bias, and flat-field corrections). Due to the 
location of SN 2009hd along a bright spiral arm of M66, contamination of the photometry 
by the host-galaxy light was a serious problem. We therefore used the template-subtraction 
technique to remove this background and hence to measure more accurately the SN magni- 
tudes. The template images of M66 were obtained with the Kitt Peak National Observatory 



13 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities for Research in Astronomy, Inc. (AURA), 
under cooperative agreement with the National Science Foundation (NSF). 
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2.1 m telescope, unde r 1" seeing conditions, about seven years prior to the SN discovery 
(IKennicutt et al.ll2003l ). The instrumental magnitudes of the SN and the reference stars in 
the SN field were measured in the subtracted i mages using th e point-spread function (PSF) 
fitting technique with the DAOPHOT package Jstetsonlll987n within IRAF. 



In order to calibrate the instrumental magnitudes to the standard photometric system, 
we used the magnitudes and colors of the local sequence stars in the SN field (Figure [1] 
and Table [T]). These values were measured from BVRI images of NGC 3627 obtained 
under photometr ic conditions on 1999 December 17 with the 1.5 m telescope at Palomar 
Observatory (see I Van Dyk et al.ll2000l ). 



SN 2009hd was observed during the first ~20 d after discovery, but it subsequently 
became too close to the Sun in th e sky. It was ob serve d again five months later (Table [2]). 
From the discovery dates given by iMonardl (120091 ) and lYamaoka fc Itagakil ( 120091 ) , we con- 
strained the explosion date to be 2009 June 19 ± 3, or JD 2,455,002 ±3. From the light 
curve s of a sample of SNe II-L, it takes ~ 9 d from explosion to B maximum (ILi et al. 



201ll ). The BVRI light curves for SN 2009hd are shown in Figure [21 relative to B maximum 
(2009 June 28 ± 3, or JD 2,455,011 ±3). Fo r comparison, we also show the light curves of 
the SNe II-L 199 0K JCappellaro et al.lll995h and 2 009kr jElias-Rosa et al.l boioh . and the 



SNe II-P 1999em (iHamuv et al.lboOll : iLeonard et alihooi lElmhamdi et al. 



20031 ) and 1992H 



( IClocchiatti et al.lll996l ). From the light-curve comparison, it appears that SN 2009hd does 
not follow a standa rd SN II-P evolution and does not behave like any SN lib (e.g., SN 1993 J; 
Barbon et al.lll995l not shown in the figure). Instead, it more closely resembles the SN II-L 
post-maximum decline. During early epochs, the light curves show a steep, almost linear, 
behavior. A linear fit to the data for the first 12 d results in decline rates f3f 2 = 0-62 ± 0.05, 
fib = 0.18 ± 0.03, /3f 2 = 0.14 ± 0.03, and (3[ 2 = 0.17 ± 0.02 mag (12 d)~ x in each band. 
We note the presence of a "bump" in the BVR light curves at ~ 10-2 d, which seems 
to follow the "shoulder" exhibited by SN 1990 K (ICappellaro et al.lll995l ). This light-curve 
feature appears in the case of several SNe II-L (IPatat et al.lll994f ). After ~ 100 d following 
B maximum, the luminosity decline settles down at a significantly lower rate. This behav- 
ior is consistent with complete trapping of the energy release from the radioactive decay 
of 56 Co to 56 Fe. In the phase range 140-260 d the decline rates are 7^ = 0.57 ± 0.08, 
7^ = 1.19 ± 0.10, and 7 J = 1.44 ± 0.12 mag (100 d) _1 , which are similar to those found for 
SN 1990K JCappellaro et al.lll995h . 



Figure [3] shows a comparison between the evolution of the absolute V magnitude 



M v , of SN 2009hd comp ared wit h the SNe II-L 1980K (IBarbon Ciatti. & Rosino 



Buta 1982; Tsvetkov 



1994aj ( IBenetti et al. 



1983), 1986E ( ICappellaro et al.lll990f ). 1990K (ICappellaro et al 



1982 



1995 



1998h . 2001cy (unpublished KAIT data), and 2009kr JElias-Rosa et al. 
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P 1992H ( 



Elmhamdi et al 



2010), wit h the SN lib 1993J (Barbon et all Il995h. and with the well-studied SNe II- 



Clocchiatti et al.l Il996h and 1999em flHamuv et all l200ll ; iLeonard et al.l 12002 



20031 ). The comp arison SNe have been corrected for extinction using pub- 



lished estimates and assuming the ICardelli. Clayton, fc Mathid ( 119891 ) extinction law. For 
SN 2009hd we have adopted a total extinction of Ay = 3.80 ± 0.14 mag (see §2.3p and 
a distance to M66 o f 9.4 Mpc (derived from the extinction-corrected Cepheid distance; 
Freedman et al.ll200ll ). Again, it is possible to see the remarkable similarity between the 
SNe II-L and SN 2009hd, primarily with SNe 1990K and 2009kr, implying that this SN also 
could be a member of this rare subclass of CC-SNe. 

Furthermore, the absolute (B — V)o, (V — R)q, (R — I)q, and (V — I)q color evolution of 
SN 2009hd are typical of a SN II-L (see Figure H]). The color curves show a rapid evolution 
from blue to red in the earlier epochs, except for the (R — J)o color curves, which initially 
exhibit a color decrease for 2-3 d, probably due to the i?-band "bump." At late times, 
however, the color becomes significantly more blue than that of other SNe II. Only the 
(R — I )o color curves for some SNe lib, such as SNe 1993J or 2008ax (not shown in the 
figure; iTaubenberger et al.l 120 111 ), have a similar behavior at such epochs. This apparent 
bluer color SN 2009hd is possibly caused by light from neighboring young stars contaminating 
our late-time measurements when the SN is faint. 



We therefore conclude that SN 2009hd displays the photometric behavior of a SN II-L. 
We caution, however, that a late plateau may have evaded detection, as a result of the absence 
of data between 30 d and 100 d. This could imply that SN 2009hd is actually a hybrid eve nt 
with characteristics of SNe II-L and SNe II-P, such as SN 1992H (jClocchiatti et al.l 119961 ). 



2.2. Spectroscopy 



Medium-resolution optical spectra of SN 2009hd were obtained on 2009 July 9 with the 
Magellan/Clay 6.5 m telescope using the Magellan Echellete (MagE; iMarshall et al.l 120081 ) 
spectrograph, and on 2010 February 15 with t he 10-m Keck-II telescope using the Deep 
Imaging Multi-Object Spectrograph (DEIMOS: LFaber et alJl2003h . 



Standard CCD processing and spectrum extraction of the MagE data were accomplishe d 
with IRAF routines. The data were extracted using the optimal algorithm of iHornd (119861 ). 
Low-order polynomial fits to calibration-lamp spectra were used to establish the wavelength 
scale. Small adjustments derived from night-sky lines in the object frames were applied. 
Fo r the MagE sp ectra, the sky was subtracted from the images using the method described 
by iKelsonl (120031 ). We employed IRAF and our own routines to flux calibrate the data 
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and remove telluric lines using the well-expose d continua of spectrophotometry standards 
f lWade fc Hor^e] Il988l : iFolev et all I2OO3L l2009h . CCD processing of t he DEIMOS spectr a 



was performed with a modified version of the DEEP pipeline (e.g., iWeiner et al.l 120051 ). 
This produced rectified, sky-subtracted t wo-dimensio nal spectra from which one-dimensional 
spectra were then extracted optimally (jHornelll986l ). The absolute flux calibration of the 
spectra was checked against the photometry and, when necessary, the spectra were rescaled. 
After that, the typical deviation of the spectra from photometry is < 10% in all bands. We 
note, however, th at the early-time spectrum was observed ~ 45° from the parallactic angle 
( jFilippenkol Il982l ) at airmass 2.1, and the flux of the late-time spectrum might be slightly 
too low because of the presence of a dead column in the CCD next to the trace of the SN; 
consequently, the shape of the continuum in both spectra could be affected. 

A near-infrar ed (NIR) spectrum was obtained with the Palomar Hale 5 m telescope 
using TripleSpec ( jHerter et al.l 120081 ) on 2009 July 15. After standard CCD processing 
and before extraction, the contribution of the night-sky lines was removed from one two- 
dimensional NIR spectral image, subtracting another two-dimensional spectral image with 
the SN placed in a different position along the slit. In the NIR, the wavelength calibration 
was performed using the OH sky lines, and a telluric absorption correction was derived from 
observat ions of the AOV star HR 6944 a nd applied to the SN data using the IDL package 
xtellcor f IVacca. Cushing. &: Raynerll2003l ). 



Figure O shows the final optical and NIR spectra of SN 2009hd. The most notable 
feature of the early optical spectrum (Figure El top panel) is the very red continuum, due to 
the appreciable extinction of this SN (see §2.3p . The spectrum also exhibits relatively weak 
spectral features, mainly due to metals, and a broad Ha emission component. The absence 
of a standard P-Cygni absorption component of Ha has also been observed in some other 
SNe II-L, such as SNe 1979C and 1980K, and this has been postula ted as a charac t eristic that 
distin guishes spectroscopically the SNe II-L from the SNe II-P ( jSchlegell Il996t iFilippenko 
1997h . 

An intriguing feature is present on the blue side of Ha, at about 6300 A. Similar absorp- 
tion i s also visible i n early-time spectra of SNe 1999em (IBaron et al.ll2000l ; iDessart &: Hillier 
20051 ) and 2005cs (IPastorello et al.l 120061 ). This feature has been interpreted either as a 
high- velocity Hi line, due possibly to an outer H-rich shell, or as Sill A6355 absorption. 
Unfortunately, the lack of spectra at early times does not allow us to definitively determine 
which of these explanations is correct. However, since the expansion velocities of the ejecta 
(~ 6400 and 3600 kms -1 , respectively) for the Ha and the putative Sill features for SN 
2009hd are similar as those for SN 2005cs at approximately the same epoch (see Figure 6 of 
Pastorello et al.ll2006l ). we surmise that the absorption near 6300 A could be more likely due 
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to Si II, rather than to high- velocity hydrogen. 



In the late-time spectrum, the continuum has become much bluer than in the early-time 
spectrum, quite probably as a result of contamination by nearby, young stars, which would 
have been included within the spectrograph slit, given the seeing conditions. The spectrum 
shows weaker Ha, [Fell] A7155, Call AA7291, 7324, and [Call] AA8498, 8542, 8662 emission 
lines (see Figure [5]). Since the SN is very close to an Hll region, which was also included in 
the slit, narrow emission lines of Ha, H/3, [O ill] A5007, [Nil] AA6548, 6584, and [S il] AA6717, 
6731 appear superposed on the spectrum. 

Finally, the spectrum in the bottom panel of Figure [5] may well be the first-ever NIR 
spectrum obtained for a SN II-L. We identified line s by comparing the spectrum to that 
of SNe II-P such as SN 2002hh JPozzo et al.l l2006h . The SN 2009hd spectrum is domi- 
nated by lines of the hydrogen Paschen series. The wavelength regions ~1. 35-1. 50 /im and 
~1. 80-2. 00 /iin were compromised, however, by strong telluric absorption. The relatively 
featureless nature of th is spectrum is reminiscent of the early-time NIR spectra of SN 1993J 
jMatthews et alJbooj ). 



The SN 2009hd optical sp ectra are shown in Figu re El relative to the spectra at simila r 
epochs of the SNe II-L 1990K JCappellaro et alJll995h and 2009kr fclias-Rosa et al.l boioh 



as we ll as of the S N lib 1993 J flBarbon et all 11995^ the SNe II-P 1992H flClocchiatti et al.l 



19961 ) and 1999em ([Leonard et al.l 120021 ) . and the SN II 200 Icy (unpublished data). All of 
the spectra are shown relative to B maximum; they have been corrected for extinction using 
published estimates and deredshifted using values from NED0. At all times, SN 2009hd 
clearly differs from SNe 1993J, 1992H, and 1999em, and more closely resembles SNe 2001cy 
and 2009kr, which also exhibit relatively weak lines (such as the weak Ha absorption trough). 
SN 1990K, also a SN II-L, shows overall stronger features than SN 2009hd. At late times, 
the features present for SN 2009hd are similar to those of normal SNe II-P, except that they 
are weaker or nearly absent. The fact that the continuum is dominated by relatively narrow 
Call and [Call] emission may imply the obscuration of the SN progenitor by dust, as well 
as the dest ruction of dust in t he circumstellar m edium (CSM); see for example the case of 
SN 2008S JPrieto et al.l l2008[ ISmith et al.l l2009h . We also compared SN 2009hd with the 
prototypical Type Il-narrow (Hn) SN1988Z (not shown in the figure); they appear similar 
blueward of Ha, but quite different on the redward side. Thus, once again we emphasize the 
resemblance of SN 2009hd to other SNe II-L. 



NED, NASA/IPAC Extragalactic Database; http://nedwww.ipac.caltech.edu/ 
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2.3. Extinction 



To derive consistent estimates of the extinction, Ay, toward SN 2009hd, we adopted 
two different methods, all based on the comparison of the SN's spectral energy distribution 
(SED) and luminosity with those of other "standard" CC-SNe. The estimates we derive 
below relate to the SN at early times. There is some hint, from both the color curves and 
spectrum, that the extinction decreased at late times; however, as discussed above, what we 
consider the more likely explanation is contamination of the SN light by neighboring stars. 

(i) We matched simultaneously the absolute (B — V) , (V — R)o, (R — I)o, and (V — 1) 
colors curves of SN 2009hd with those of other SNe II-L (see Figure H]). The best match is 
for Ay = 3.9 ± 0.1 mag. This comparison was primarily weighted by the larger number of 
data points at early epochs. 

(ii) The observed extinction was estimated by comparing the early-time optical SED of 
SN 2009hd with those of SNe 1990K, 2001cy, and 2009kr at similar epochs. Spectra of the 
comparison SNe were first corrected for redshift and extinction, and scaled to the distance 
of SN 2009hd. The average of the good matches in all cases is Ay = 3.7 ± 0.2 mag. 

We find that these extinction estimates from the two quite similar methods agree rea- 
sonably well with each other. Hence, we adopt Ay* = 3.80 ± 0.14 mag as the extinction 
toward SN 2009hd. The foreground Galactic c omponent to this extinction is a relatively low 
0.11 mag (jSchlegel. Finkbeiner. fc Davislll998l ). so, therefore, most of the extinction must be 
local to the SN. 



This extinction estimate is consistent with the constraint on the reddening we can 
derive from the Magellan+MagE spectrum. We measured the equivalent width (EW) of the 
Nal D line at the host-galaxy redshift (zo — 0.0024) and found for the two line components, 
EW(Nal Dl A5890) = 1.25 ±0.05 A and EW(Nal D2 A5896) = 1.10 ±0.05 A. Considering 
that EW(Nal Dl) hr + EW(Nal D2) hr = EW(Nal D) lr (Elias-Rosa et al., in prep.), where 
"hr" and "lr" denote "high" and "low resolution," respectively, this results in EW(Nal D) 
= 2.35 ± 0.07. However, since the rat io between th e Nal D component lines is close to 
saturation (1.14, with saturation at 1.1; ISpitzerl Il948h . and given the con siderable scatter in 



th e distribution of SNe i n the EW(Nal D) vs. E(B - V) plane (e.g., see felias et aDl200# 5 



or 



Poznanski et al.ll201ll ) and the large variance in this relation found in the literature, the 
estimated value of E(B — V) tot can only be considered as a lower limit. For example , using 



the relation between reddening and EW(Nal D) (for all type of SNe) from lElias et al.l ( 120071 . 
E(B - V) = 0.11 x EW(Nal D) -0.01), and assuming the Cardelli et al. (1989) reddening 



15 http://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=60. 
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l aw, with updated bandpass wavelen gths, and a Galactic foreground E(B — V) = 0.03 mag 
JSchlegel. Finkbeiner. fc Davislll998h . we estimate E(B - V) tot £ 0.29 ± 0.03 mag (A v £ 
0.9 ± 0.09 mag and E\V — 7] tot ^ 0.35 ± 0.06 mag). If instead, we use the relation given 
bv loiivares E. et all feoioh for SNe II-P (E(B - V) « 0.25 x EW(Nal D)), we estimate 
i?(.B — V) tot ^ 0.62 mag. In both cases, these lower limits are consistent with the two direct 
extinction estimates above. 

We note here the possible temporal evolution of the Na I feature. While it is strong at 
early times, indicative of the high extinction toward the SN, the feature's strength appears 
to be reduced by ~ 20% (EW[Nai D] = 1.87 ± 0.05) at late times. Such variability in the 
Nal D lines has been observed before — for examp le, in the faint transient SN 2008S in 



NGC 69 46 (ISmith et al 



2009 



Botticella et al.l 120091 ). or in thermo nuclear SNe, such a s SN 



2006dd f lStritzinger et aUEOllt ). and SN 1999cl and SN 2006X (see lBlondin et al.l 120091 . and 
references therein). However, in these latter two cases, the EW of the lines increased with 
time, rather than decreased. Considering the expected recombination of Na with time, such 
a variation for SN 2009hd, if real, may be due to ionization changes in the CSM, as a result 
of the weaker radiation field from the SN at later times. This effect could also be related to 
the possible decline in the local extinction mentioned in §2.21 On the other hand, another, 
possibly more likely, explanation for the apparent change in the Nal D EW is dilution of 
this spectral feature by light from nearby stars and H II regions, as previously noted. 



3. Identification of the Progenitor Candidate 



Shortly after the discovery of SN 2009hd, we isolated 12 archival HST/WFPC2 CR- 
SPLIT image pairs of M66 in F555W (~ V) and 5 pairs in F814W (~ J), in order to identify 
any stars in the images that could be candidates for the SN progenitor. These images were 
obtained between 1997 November and 1998 January by program G O-6549 (PI: A. Sa ndage) 



and used originally to derive a Cepheid-based distance to M66 (ISaha et al.l Il999f ). The 
exposure times in each band were 2400-2500 s. Combining the images in each band over 
all epochs of the Cepheid observations yields a total image depth of ~ 16 hr in F555W and 
~ 7 hr in F814W. The combination of the pre-explosion images for each band scaled t o 



the PC chip pixel scale (0'/045 pix -1 ), using the drizzle algorithm (IFruchter fc Hookll2002l ). 
is available from the HST Legacy Archive 17 !. The input images in the drizzle algorithm 
are weighted according to the statistical significance of each pixel, and drizzle removes the 



16 http: / / archive.stsci.edu/hst/. 
17 http://archive. stsci.edu/hlsp/. 
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effects of geometric distortion on both image shape and photometry. In addition, it combines 
dithered images in the presence of cosmic rays and improves the resolution of the mosaic. 
Thus, the use of the drizzled image mosaics for the identification of a progenitor candidate 
is ideal. 

Photometry of the post-explosion HST images was performed using the DOLPHOTr^l 
package especially designed for ACS. The post-explosion images were all obtained at the 
same epoch in each band and with a small dithered offset between exposures. Thus, we 
measured these relative offsets with respect to one fiducial image before running DOLPHOT. 
The output from the package automatically includes the tr ansformation from flight-system 



F555W and F814W to the correspond ing Johnson- Cousins (IBesselll I199CH ) magnitudes (in V 



and /), following the prescriptions of ISirianni et al.l (120051 ) . We include photometry of the 
SN in Table H 

For the SN progenitor candidate identification, it was necessary that we use high- 
resolution images of the SN, ideally at the same or better resolution as the pre-SN HST 
image. Accordingly, we obtained HST Advanced Camera for Surveys (ACS; Wide Field 
Channel, ~ 0'.'05 pix -1 ) images in F555W and F814W on 2009 December 14, as part of our 
Target-of-Opportunity program GO- 11 575 (PI: S. Van Dyk). The SN was still quite bright 
at the time of these observations (see Table [2]), although it was not saturated in the images. 
The observations in each band were in pairs of dithered 260 s exposures. These individual 
exposures were drizzled to produce final mosaics in each band. A comparison of a HST SN 
image to a pre-SN image is shown in Figure [7J 

We achieved high-precision relative astrometry between the pre-explosion drizzled WF2 
image in each band (the SN site is located on this chip of the pre-explosion images; this 
image is of the same depth as the drizzled full mosaic image) and the post-explosion ones by 
geometrically transforming the former to match precisely the latter. Using the IRAF task 
geomap, we carried out several geometric transformations between the two sets of images, 
using different numbers of point-like sources in common between the two datasets (from ~10- 
30 manually selected sources to ^ 300 sources using an automated matching algorithm), all 
with root-mean-square (rms) uncertainty <^ 0'.'15. We then transformed the pre-explosion 
images to the post-explosion ones using the IRAF task geotran. The positions (and their 
uncertainties) of the SN and the progenitor candidate were generally obtained by averaging 
the measurements from two centroiding methods, the task daofind within IRAF/DAOPHOT 
and imexamine within IRAF. For all solutions, the resulting SN position in the pre-explosion 



IS 



The ACS module of DOLPHOT is an adaptation of the photometry package HSTphot. We used vl.l, 
updated 2010 January 6, from |http : / /purcell. as . arizona.edu / dolphot /\ 
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images in the two bands were all consistent. 

As a result, in the pre-SN F814W image we identify an object very near the SN position 
which we consider to be the progenitor candidate. It is not detected at this position in the 
F555W image; see Figure [HJ The average final pixel position in the pre-SN F814W image 
is [1074.82, 1634.13] for the SN, and [1074.43, 1634.10] for the candidate in the transformed 
SN image. The candida te is o n e of t he detected stars near the edge of the pre-explosion 



error circle suggested by iMaozl (120091 ) in his initial attempt to identify the SN progenitor. 
The differences between the SN and the progenitor candidate positions, compared with the 
total estimated uncertainty in the astrometry, are given in Table |3j This latter uncertainty 
was calculated as a quadrature sum of the uncertainties in the SN and progenitor candidate 
positions, and the rms uncertainty in the geometric transformation. 

From the results in Table [31 it can be seen that the difference between the SN position 
and the position of the progenitor candidate in F814W is larger in right ascension than 
the measurement uncertainties. This offset could be due to the proximity of one or more 
sources to the progenitor candidate position, which may have affected the measurement of 
the object's centroid, or to the difficulty in making accurate centroid measurements, given 
the complex background in the bright spiral arm of the host galaxy on which the progenitor 
candidate is located. In Figure |8]one can see that the reference object, labelled "B" (~ 0'.'25 
northwest of the candidate), is visible in the F814W image, but is essentially invisible in 
the F555W image. The opposite is true for the object, labelled "C" (~ 0'.'33 southwest 
of the candidate). We emphasize that no other point-like source can be isolated within a 
1.5cr (i.e., <2 pixel) radius from the progenitor candidate position (see Figures [7] and |8]) . 
The alternative possibility is that we have not detected the progenitor for SN 2009hd. We 
consider both of these possibilities below. 



4. The Nature of the Progenitor 

Next, we attempted to measure the brightness of the progenitor candidate. We per- 
formed photometry on each of the individual pre-explosion HST images in both bands using 
HSTphot0. (HSTphot is not designed to run on the drizzled image mosaics.) We ran the 
package on the CR-SPLIT pairs for both bands, with option flag 10, which combines turn- 
ing on the local sky determination, turning off empirically determined aperture corrections 
(using default values instead), and turning on PSF residual determination, with a total 



19 HSTphot is a stellar photometry package designed for use with WFPC2 images (jDolphirJ 120001 ) . We 
used vl.l.7b, updated 2008 July 19. 
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detection threshold of 3a. As a result, we obtained mp8i4w = 23.54 ± 0.14 mag for the 
progenitor candidate. HSTphot did not detect in F555W the candidate detected in F814W, 
even with the detection threshold set as low as la. To be sure that this object is truly not 
detected in F555W, we repeated the measurements, varying the number and order of the 
images for each filter input into HSTphot (we note that the images in each band can be 
split into two subsamples, depending on the various camera pointings), both automatically 
and forcing HSTphot to make a detection at the centroid position of the object in F814W. 
Complicating the detection is the background. The number of counts around the SN position 
in the F555W pre-SN images is nearly a factor of two higher than in a similarly sized region 
measured in a "empty" area, relatively free of galaxy contamination. 

We note that HSTphot reported the progenitor candidate as having an "object type" flag 
of "1," indicating that the source is likely stellar. Additionally, the object had a "sharpness" 



point- like: 


Dolp 


lin 


2000) 


(Leonard et al. 


2008). E 



1.5; hence, it is reasonable to call the 



candidate a "good and clean" star (x < 1-5, see iDolphinl 120001 ) . We therefore consider it 
very likely that the candidate is a single star. 

Of course, given the offset of the SN position from this star, we must still take into 
account the possibility that we have not actually detected the progenitor in either band. 
We therefore attempted to constrain the brightness of an undetected progenitor by adding 
an artificial star at the precise SN position (see § E]) in the pre-SN images for both bands. 
The PSF for the input artificial star was created from bright, isolated stars seen in the 
drizzled image mosaic in each band using IRAF/DAOPHOT. The brightness of the artificial 
star was varied, until it was practically not detectable above the background. Then, the 
instrumental magnitudes of the artificial star were measured in the images using PSF fitting, 
and subsequently calibrated to the flight-system photometric system using the magnitude 
and colors of 8 brighter stars in the vicinity of the SN, measure with HSTphot in the HST 
WFPC2 images. We find that > 26.1 mag in F555W and > 23.5 mag in F814W. Note that 
the limit in F814W is very similar to the measure brightness of the progenitor candidate in 
this band. Correcting by the same total extinction as assumed for the SN ( Ay* = 3.80±0.14 



mag: see § 12. 3p and the adopted distance modulus (fi = 29.86 ± 0.08 mag; iFreedman et al. 



20011 ) . we find a limit on the absolute magnitude of an undetected progenitor, M v ^ —7.6 



mag. Furth ermore, assum ing a bolometric correction of —0.29 mag for supergiants at solar 



metallicity ( iKuruczl 119931 ; Kurucz Atlas 9 models, CD-ROMs 13, 18; see below for the 
metallicity measurement in the SN environment) appropriate for (V — I)q ^ 0.99 mag (this 
is the limit on the color that arises from the reddening-corrected upper limits at V and /, 
above), we find an upper limit to such a star's bolometric luminosity of Lboi ^ 10 5 04 L . 
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The limit on the effective temperature, from the color limit, is then T e g <; 5200 K. We place 
this limit on the progenitor in the Hertzsprung- Russell (HR) diagram in Figure [91 

We can compare the limit shown in the figure t o massive-star evolutionary model s 



without rotation and with rotation (t> ini = 300 kms -1 iHirschi. Meynet. fc Maeder 1 12004 ). 
One can see that the appropriate initial masses to consider, of those that are available for 
these models, are 15 and 20 M . (We note that the majority of stars do not rotate as 
fast as those in the models used here; e.g., less tha n about 10% of yo ung Galactic O-type 



stars have rotation speeds of 300 km s _1 or more; iPenny et al.ll2004 ) We also show the 
evol utionary tracks for a mo del star of initial mass 20 M with pulsationally enhanced mass 
loss (jYoon &: Cantielloll2010l ). We assume that these tracks correspond to solar metallicity. 
We determined directly the metallicity in the SN environment from a neighboring H II region 
detected in the late-time SN spectrum (see § I2.2p . Measuring the Ha and [N II] A6584 line 
intensities and app lying a cubic-fit relati on between the ratio of these two lines and the 
oxygen abundance (jPettini fc Pagell 12004 ). we find 12 + log(0/H) = 8.43 ± 0.05, which is 
somewhat lower t han the solar abunda nce, 8.69 ± 0.05 (jAsplund et al.ll2009i ). We adopt the 
criterion given by ISmartt et al.l (120091 ). which, taking into account the uncertainties in the 
different metallicity estimates, considers 12 + log(0/H) = 8.4 as the dividing line between 
solar and subsolar metallicity evolutionary tracks. 

To convert this progenitor detection limit to an upper limit on the initial mass of the 
progenitor, we considered the luminosities of the evolutionary tracks in the post-red super- 
giant (RSG) phase (see Figure [9]). We made this estimation three different ways. First, we 
compared the lum inosity upper limit to the stellar endpoints of the evolutionary tracks, as 
Fraser et all J2009h did in the case of the SN II-L 2009kr. Doing so, we find that M ini £ 2OM 
for the models with no rotation, and <; 15 M & for the models with rotation. If we, instead, 
consider where the upper limit intersects the evolutionary tracks, we find a similar upper 
initial mass limit, <^ 20M Q , for the non-rotating models, and <^ 18 M for the models 
with rotation. If we consider the model with pulsationally-induced mass l oss at 20 M^, one 
can se e that it naturally achieves terminus blueward of the RSG phase. lYoon fc Cantiello 
(120101 ) found that their 17 M model only experiences pulsations after Ne burning. One can 
imagine, then, that, as the models approach 20 M from that lower mass, that pulsations 
become a stronger effect, and the models also terminate toward the blue, rather than the red. 
Thus, conservatively, we can say that, if pulsations are important for these massive stars, 
the luminosity limit could correspond to a mass limit <^ 20 M ; however, more liberally, 
this can extend to somewhat a lower mass limit of <^ 18-19 M & . 
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Discussion and Conclusions 



Observations at early times and at later (nebular) phases suggest t hat SN 2009hd is 



a SN II-L, similar to SN 1990K flCappellaro et al.lll995l ) and SN 2009kr flElias-Rosa et al. 
2010l ) . However, the lack of data between ~ 30 and 100 d leaves open the possibility that 



SN 2009hd may be a hybrid between SNe II-L and SNe II-P. If it is a SN II-L, it is one of 
the most extinguished examples detected thus far, possibly as a result of proximity to a dust 
lane that is clearly evident in pre-SN HST image s. Projection effects may further play a 
role, given the inclination of M66 (65°; lTullylll988l ). 



Comparing HST images of the host galaxy prior to explosion and of the SN at late times, 
the resulting high-precision relative astrometry has provided us with a means of identifying 
the SN progenitor. A progenitor candidate is possibly identified in the available F814W, 
but definitely not in the F555W, pre-SN images. Even with the high spatial resolution of 
the HST WFPC2 images, the uncertainties in the astrometry are large enough that we are 
unable to make a definitive identification. Given the relative proximity of the host galaxy 
to us (9.4 Mpc), this is rather surprising. The high extinction, combined with the complex 
background in the SN environment, likely lead to this level of uncertainty. 

We can constrain the properties (luminosity, effective temperature) of the progenitor 
star for SN 2009hd. Although the inferred properties are certainly consistent with a red 
supergiant (RSG), the constraints also allow for the possibility that the star is more yellow 
than red. A yellow color could imply that it exploded in a post-RSG evolutionary state (or 
perhaps that it is an RSG in a binary system with a bluish companion). The progenitor may 



have evolved back to the blue a: 



mass (e.g.. iPodsiadlowski 



1992 



ter reaching the RS G stage, for example, because it has lost 



Mevnet et al .1120 111) as a result of rotationally induced winds 



Hirschi. Meynet. fc Maeder 1 120041 ). or through pulsationally enhanced mass loss (e.g., 



Yoon fc Cantielldl2010l ). The fact that the H lines are weak in the nebular spectrum may 
indicate that the star exploded with a low-mass H envelope, presumably lost through high 
mass loss via stellar winds prior to the termination of He burning. Conseq uently, the star 



would evolve toward the blue before its terminus (e.g., iMeynet et al.ll201lL and references 



therein). Through comparison with high-mass stellar evolutionary tracks at solar metallicity, 
we are able, conservatively, to constrain the initial mass of the progenitor to Mi n ; <^ 20 M . 

This limit to the initial mass estimate for the SN 2009hd proge nitor is consistent with 
the ran ge in mass found f or the SN II-L 2009kr, M ini « 18-24 M felias-Rosa et al.lboiO 



see also iFraser et all 120091 ) . The progenitor of SN 2009kr was inferred likely to be a yellow 
supergiant. The progenitor of the unusual, high-lumin osity SN II-P 2008cn a lso appears to 
have been a yellow supergiant, with M ini = 15 ±2 M Q felias-Rosa et al.ll2009l ). SNe 2009hd, 
2009kr, and 2008cn may be a bridge between the RSGs that explode as SNe II-P and the 
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high-mass, lumino us blue variables or more extreme RSGs that may explode as SNe Iln 
jSmith et al.lboilh . 

Ultimately, a very-late-time set of multi-band images of the SN 2009hd field should be 
obtained with HST and the Wide Field Camera 3 (WFC3) many years later, when th e SN 



has greatly faded, as has been done for the SN Hn 2005gl (|Ga. 



for the SN lib 1993J and SN II-P 2003gd flMaund fc Smarttl 



-Yam fc Leonard! 120091 ) and 



20091 ). to better distinguish 



the constituents of the SN 2009hd environment and help decipher the true nature of its 
progenitor. 
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Fig. I.— /-band image of SN 2009hd in NGC 3627 (M66) obtained with the 1.3 m SMARTS 
telescope+ANDICAM at CTIO on 2009 July 10 (field of view ~ 6' x 6'). The sites of 
SNe 1973R, 1989B, and 1997bs, as well as the local photometric sequence stars, are indicated. 
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Fig. 2. — Le/f: Early-time ( ^ 40 d after 5-band maximum) optical light curves of SN 2009hd, 
together with those of the SNe II-L 1990K (dot-dashed lines) and 2009kr (solid lines), and 
of the SNe II-P 1999em (dashed lines) and 1992H (dotted lines). Right: Optical light curves 
of SN 2009hd, including late-time data. The comparison light curves are adjusted in time 
and magnitude to match those of SN 2009hd. 
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Fig. 3. — Absolute V light curve of SN 2009hd (filled circles) along with those of the SNe II- 
L 1980K (empty squares), 1986E (empty stars), 1990K (5-pointed stars), 2001cy (dashed 
lines), and 2009kr (solid lines); the peculiar SN II-L 1994aj (3-pointed stars); the SN lib 
1993 J (dotted lines); and the SNe II-P 1999em (empty circles) and 1992H (empty triangles). 
Distances and extinction estimates for all SNe, except SN 2009hd, have been adopted from 
the literature. The magnitudes of SN 2009hd have been corrected for the assumed red- 
dening, E{B — V) = 1.23 mag, and the adopted distance modulus, fi = 29.86 ± 0.08 mag 
( IFreedman et al.ll200ll ). 
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Fig. 4.— Intrinsic (B - V) , (V - R) , (R - J) , and {V - f) color curves of SN 2009hd, 
compared with the intrinsic color evolution of the SNe II-L 1990K (5-pointed stars), 2001cy 
(dashed lines), and 2009kr (solid lines); the peculiar SN II-L 1994aj (3-pointed stars); the 
SN lib 1993 J (dotted lines); and the SNe II-P 1999em (empty circles) and 1992H (empty 
triangles). Distances and extinction estimates for all SNe, except SN 2009hd, have been 
adopted from the literature. The color of SN 2009hd has been corrected for the assumed 
reddening, E{B — V) = 1.23 mag, and the adopted distance modulus, /i = 29.86 ± 0.08 mag 
( IFreedman et al.ll200ll ). 
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Fig. 5. — (a) Optical spectral evolution of SN 2009hd. The late-time spectrum appears 
to have been contaminated by light from stars near the SN. (b) A near-infrared spectrum 
obtained using TripleSpec with the Palomar 5 m Hale telescope on 2009 July 15. Spectra 
in both panels have been corrected for their host-galaxy recession velocities, but not for 
extinction. Ages are relative to B maximum light. The locations of the most prominent 
spectral features are marked. 
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Fig. 6. — Early and late-time optical spectra of SN 2009hd (the top spectrum in each 
panel, shown in bold), along with comparison spectra of the SNe II-L 1990K and 2009kr, 
SN lib 1993J, SNe II-P 1992H and 1999em, and SN II 2001cy. All spectra have been corrected 
for their host-galaxy recession velocities and for extinction (see text). Ages are relative to 
B maximum light. 
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Fig. 7.— (a) Subsections of the pre-explosion HST WFPC2/F814W images of M66, and (b) 
the post-explosion HST ACS/F814W image of SN 2009hd. The positions of the candidate 
progenitor and the SN are indicated by tick marks. 
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Fig. 8. — Subsections of the pre-explosion HST WFPC2 images of M66 in (a) F555W and (b) 
F814W. The position of the SN is indicated by a 5a positional uncertainty ellipse (0'.'04 x 0'.'03 
for the F555W image and 0'.'08 x 0'/04 for the F814W image). The four objects "A," "B," 
"C," and "D" were not used in the image astrometry, due to their relative faintness; they 
are labeled here only to aid the eye in locating the progenitor in the field. 
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Fig. 9. — A Hertzsprung- Russell diagram (T e g (K) vs. Lboi (L©)) showing the upper limit 
considering a contaminated progenitor (filled 5-pointed star). Model stellar evolutionary 
tracks for solar metallicity are also illustrated for M ini of 15 and 20 M without rotation 
(short- dashed, and long-dashed gray line s) and with rotation of = 300 k ms -1 (dot-short- 
dashed, and dot-long- dashed black lines; iHirschi. Meynet. fc Maeder I I2004T) : and for 20 M© 



consi dering pulsation-induced mass loss enhancement (solid black line; 
2010h . 



Yoon fc Cantiello 
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Table 1. Adopted Magnitudes for the Local Sequence Stars in Figured] 



Star 


B 


V 


R 


I 




(mag) 


(mag) 


(mag) 


(mag) 


1 


13.72±0.01 


13.17±0.01 


12.83±0.01 


12.50±0.01 


2 


16.26±0.01 


15.47±0.01 


14.97±0.01 


14.49±0.03 


3 


17.16±0.01 


16.54±0.01 


16.16±0.02 


15.77±0.01 


4 


16.77±0.01 


16.24±0.01 


15.92±0.01 


15.55±0.01 


5 


16.56±0.01 


15.95±0.01 


15.59±0.01 


15.18±0.01 


6 


18.03±0.01 


16.49±0.01 


15.41±0.01 


13.49±0.24 
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Table 2. Photometry of SN 2009hd 



Date 


JD 

-2,400,000.00 


B 

(mag) 




V 
mag) 






R 

^mag) 






I 

mag) 




Source 


06/07/09 


55019.48 


18.25±0.05 


16 


73±0 


06 


15 


82±0 


05 


14 


66±0 


04 


CTIO 


r\ o 1 f\'~7 1 f\C\ 

08/07/09 


55021.46 


18.43±0.04 


16 


78 ±0 


03 


15 


83±0 


04 


14 


69±0 


04 


CTIO 


09/07/09 


55021.70 


18.44±0.06 


16 


75±0 


03 


15 


75±0 


04 








P60 


10/07/09 


55022.70 


18.56±0.05 


16 


78±0 


04 


15 


81±0 


02 


14 


68±0 


04 


P60 


10/07/09 


55023.45 


18.41±0.07 


16 


75±0 


02 


15 


76±0 


03 


14 


71±0 


03 


CTIO 


12/07/09 


55025.46 


18.57±0.08 


16 


83±0 


06 


15 


79±0 


06 


14 


69±0 


04 


CTIO 


14/07/09 


55026.69 


18.65±0.03 








15 


82±0 


06 








P60 


16/07/09 


55029.46 


18.73±0.07 


16 


84±0 


02 


15 


83±0 


05 


14 


83±0 


03 


CTIO 


18/07/09 


55031.46 


18.98±0.06 


16 


94±0 


02 


15 


94±0 


05 


14 


81±0 


03 


CTIO 


24/11/09 


55160.71 




20 


64±0 


05 


19 


63±0 


02 


18 


75±0 


02 


LT 


14/12/09 


55179.58 




20 


87±0 


01 








18 


95±0 


01 


HST 


19/12/09 


55184.86 




21 


09±0 


10 


20 


11±0 


08 








Mag 


20/01/10 


55216.55 




20 


85±0 


04 


20 


30±0 


05 


19 


54±0 


04 


LT 


15/02/10 


55242.93 










20 


83±0 


05 


19 


97±0 


07 


Keck 


19/03/10 


55274.54 




21 


43±0 


09 


20 


99±0 


08 


20 


33±0 


13 


NOT 



Note. — CTIO = CTIO 1.3 m SMARTS telescope + ANDICAM; P60 = Palomar 1.5 m 
telescope + CCD camera; LT = 2.0 m Liverpool telescope + RATCam; HST = Hubble Space 
Telescope + ACS camera; Mag = 6.5 m Baade Magellan Telescope + IMACS; Keck = 10 m 
Keck-II telescope + DEIMOS; NOT = 2.6 m Nordic Optical Telescope + ALFOSC. 
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Table 3. SN and Progenitor Candidate Position Comparison 



Angular Quantity 


V (a/5) 


/ (a/5) 


Uncertainty in the progenitor position (mas) 




14/5 


Uncertainty in the SN position (mas) 


1/3 


1/3 


Geometric transformation (mas) 


8/5 


7/6 


Total uncertainty (mas) 


8/6 


15/8 


Difference in position (mas) 




20/1 



Note. - Uncertainties in the SN and candidate position, in mil- 
liarcsec (mas), were estimated as the standard deviation of the aver- 
age. Geometric transformation errors are derived from the positional 
differences between the fiducial stars used in the transformation. The 
total uncertainty is the quadrature sum of these uncertainties. The 
last line lists the residual difference between the SN and progenitor 
position after the geometric transformation. 



